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Abstract 

This paper outlines a method for the stress analysis of bonded composite corrugated panel facesheet 
to flange joints. The method relies on the existing HyperSizer Joints software, which analyzes the bonded 
joint, along with a beam analogy model that provides the necessary boundary loading conditions to the 
joint analysis. The method is capable of predicting the full multiaxial stress and strain fields within the 
flange to facesheet joint and thus can determine ply-level margins and evaluate delamination. Results 
comparing the method to NASTRAN finite element model stress fields are provided illustrating the 
accuracy of the method. 


I. Introduction 

Preliminary design of corrugated panels requires the rapid consideration of a wide range of panel 
configurations (materials, layups, facesheet thicknesses, stiffener geometries, etc.) to determine the 
optimum design for a given set of loads. To assess the efficiency of a given panel design, ply-level 
stresses are used in one or more failure criteria to estimate the panel design’s minimum margin with 
respect to failure, given panel level loads. The standard method of determining ply-level stresses, finite 
element analysis, however, is not well suited for the rapid consideration of (sometimes) thousands of 
potential panel designs in which the geometry is not fixed. In contrast, the HyperSizer Structural Sizing 
Software (ref. 1) uses rapid, closed-form solution methods on the level of the panel, laminate, and ply to 
evaluate panel margins and determine the optimum lightweight design. This software still relies on a 
global finite element model to determine the panel level loads, but it calculates the internal stresses within 
the panel based on efficient analytical solutions. 

The present paper is concerned with the local stresses that arise within the bonded joint that exists 
between the facesheet and stiffener flanges in a composite corrugated stiffened panel. Existing 
HyperSizer capabilities analyzed corrugated panels using the ply-level stresses in each component of such 
a stiffened panel (e.g., stresses in the flange, web, crown, and facesheet), but did not consider stress 
concentrations due to the facesheet-stiffener bond. To address this limitation, the methodology described 
herein combines an existing bonded joint analysis capability within the software with a beam analogy 
analysis of the corrugated panel geometry. The beam analogy analysis provides the boundary loads, 
which vary based on the panel configuration, to the joint analysis. The joint analysis then calculates ply- 
level in-plane and intralaminar stresses and adhesive stresses that arise in the joint between the facesheet 
and stiffeners, which can be used to determine ply-level margins. This new capability has been 
incorporated within HyperSizer, enabling sizing based on global load cases down to the level of the stress 
concentrations arising due to the bond of each stiffener in each panel. Results are presented that compare 
local stress concentrations predicted by HyperSizer to those predicted using continuum finite element 
models to represent the facesheet to stiffener flange joint of several corrugated panel configurations. 
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II. Joint Analysis Model 

The joint analysis model employed is based on Mortensen's unified stress analysis method (refs. 2 and 
3), which has been extended and implemented in the HyperSizer Joints composite bonded joint analysis 
software. Extensions include accommodation of transverse in-plane straining, hygrothermal loads, 
computation of the local in-plane and intralaminar stresses throughout the adherends, accommodation of 
pressure loading, and delamination prediction (refs. 4 to 9). Compared to other analytical (i.e., non-FEA) 
methods used for bonded joint analysis, the HyperSizer method is capable of handling more general 
situations, including various joint geometries, asymmetric and unbalanced laminates, and more general 
loading and boundary conditions. A wide range of joint types may be considered, and the adherends, 
which were originally modeled as classical laminates in cylindrical bending, are now considered to 
undergo ‘generalized cylindrical bending’, in which transverse straining is accommodated. Both linear 
and nonlinear behavior of an adhesive layer is admitted in the analysis. For linear analysis, the adhesive 
layer is modeled via a traction-separation law that responds similarly to the linear behavior cohesive 
elements within the ABAQUS finite element software (ref. 10). This can be used to simulate a physical 
adhesive layer present in the joint, or, for cases without a non-negligible adhesive layer, the traction- 
separation model represents a zero-thickness cohesive layer that can be given a high penalty stiffness to 
hold the adherends together (ref. 11). Inclusion of nonlinear adhesive behavior in the analysis is 
accomplished through the use of a secant modulus approach for the nonlinear tensile stress-strain 
relationship in conjunction with a yield criterion. 

The equilibrium equations for each joint type are derived through direct imposition of force and 
moment equilibrium on joint elements, and by combination of the aforementioned equations and 
relations, a set of governing ordinary differential equations is obtained. The governing system of 
equations is solved numerically using Mortensen and Thomsen’s (ref. 3) ‘multi-segment method of 
integration,’ yielding laminate-level fields and adhesive stresses that vary along the joint in each 
adherend. After the governing equations are solved, the ply-level in-plane stress components in the 
adherends can be calculated from Classical Lamination Theory (CLT). After solving for the in-plane 
stresses, the intralaminar stress components in the adherends are obtained through integration of the 
point-wise equilibrium equations. The details of HyperSizer’s stress analysis method are described by 
Zhang et al. (refs. 5 and 8). It is important to recognize that the HyperSizer joint analysis method is very 
efficient; the execution time for a typical problem is approximately 1/40 sec. Furthermore, problem set up 
and post-processing is straightforward, enabling rapid consideration of a wide range of joint 
configurations for sizing optimization. 

In order to model the flange to facesheet joint within a composite corrugated panel, appropriate 
boundary loading conditions must be applied within the HyperSizer Joints analysis. As shown in figure 1, 
the analysis domain for a hat stiffened panel is the flange and facesheet region, where the closed span and 
web attachments to the analysis domain are accounted for through the boundary loading conditions. This 
region is, in essence, a bonded doubler, which can be analyzed by HyperSizer Joints. A two sheet 
corrugated panel (which has continuous flanges between the hats) is shown in figure 2. The analysis 
region is again a bonded doubler. The next section describes the determination of the appropriate 
boundary loads to apply (BCs #2 and #3 in figure 1 and BCs #1 and #4 in fig. 2) within the HyperSizer 
Joints analysis. 


III. General Solution for Boundary Loads 

The goal of this section is to determine the internal forces and moments in the components of a 
corrugated stiffened panel subjected to an applied shear force (V), an applied distributed load (force per 
unit length) due to internal pressure (P), and either an applied normal force (F) and moment (M) or a 
symmetry condition resulting in a normal reaction force ( F) and moment (M) at point H, as shown in 
figure 3. Note that point H is the midpoint between the corrugations (see fig. 3). For a two sheet 
corrugated panel, which has no open span, L os is set to zero. A beam analogy model is used to determine 
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the internal forces and moments, the knowledge of which provides the unknown applied boundary 
conditions indicated in figures 1 and 2. The two separate boundary conditions (applied force and moment 
versus symmetry) are discussed in Sections III. A and III.B. It should be noted that in the following beam 
analogy model, the direction along the facesheet is referred to as the x-direction (fig. 3), while in the joint 
analysis, this is the y-direction (Figs. 1 and 2). 

In figure 3, member BD is referred to as the closed span (cs), member AB is referred to as the web 
(w), member AC is referred to as the crown (cr), member GB is referred to as the facesheet-flange 
combination (com), and member HG is referred to as the open span (os). Symmetry conditions are applied 
at points D and C that allow z-direction translation (w), but no x-direction translation (u) and no rotation. 
The model is restrained against w displacement at point B, the intersection between the hat flange and 
web. This choice was motivated by comparison to finite element analysis results for acreage deformation 
of a hat stiffened panel model subjected to pressure loading. Note that neither horizontal 
(x-direction) translation nor rotation is restrained at point B. 


BC#1 

u = v = 0i = 0 
Applied Q y 


BC#2 

w=0 

Applied N y 
Applied A/ xy , closed span 



Figure 1 . — Analysis geometry for a composite hat stiffened panel facesheet to flange joint. 


BC#1 

Applied N y , Q y 
Applied A/ xy , closed span 


BC#2 (symmetry) Applied M y , closed span 







Figure 2. — Analysis geometry for a composite two sheet stiffened panel facesheet to flange joint. 
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A. Applied Force and Moment Boundary Conditions between Corrugations 

A free body diagram of the analysis geometry is considered, as shown in figure 4. Applying global 
equilibrium, 


o 

II 

cC 

IXI 

— » b 7 +-p-v = 0 

2 

(i) 

M 

II 

o 

— ^ D x — F + C x = 0 

(2) 


2X<. 


D 


= 0 


M D +M r +hC r +-V P-M-B ,L„ =0 


l C 


8 


( 3 ) 




Figure 3. — Analysis geometry for a hat stiffened panel 
B z 
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Substituting for B z in equation (3) using equation (1), 


M D +M c + hC x + V 




+ P 


f S_ _5_ 2 ^ 


— M = 0 


( 4 ) 


Next a free body diagram of member HGBD is considered (as shown in fig. 5), and equilibrium is 
applied, 


I>z = 0 B z -BA z+ ^P-V = 0 

-> D X -F + BA X = 0 


Combining equation (5) and (1) gives, 


BA Z = 0 


and combining equation (6) and (2) gives, 


BA X = C X 


( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 


Therefore, equilibrium has established the relationship among the forces in the members, as well as one 
equation (4) relating the moments. Considering C x , M c , and M D to be the unknown quantities in the 
formulation, two additional equations are needed in addition to equation (4). These equations will be 
obtained by applying beam theory deformation equations while imposing continuity of rotations and 
x-direction displacements at point B. 

It is noted from the boundary condition at point D depicted in figure 3 that the x-direction 
displacement and the rotation of member BD at point D are zero. Considering the free body diagram of a 
section of member BD as shown in figure 6, and introducing a new x coordinate direction, 


B z 


M 

. c 

H 

q m b > 

i b ^. Ax 

D 

F C 

II 1 1 1 \y 

' i i i m ii m 


d\ Md 

r7~ D 


ba 3 


Figure 5. — Free body diagram of member HGBD (closed span, facesheet-flange 
combination, and open span). 




M cs (x) 


d^ d 

V Dx 


V cs (x) P 

Figure 6. — Free body diagram of a section of member BD. 


NAS A/TM— 2008-2 1 543 8 


5 


2> z =o 


( 9 ) 


V cs~ P X 




M d -M cs+ V cs x-P— = 0 


( 10 ) 


Combining equations (9) and (10), 


^cs 


— M / ) + P — 

u 2 


( 11 ) 


Beam theory (ref. 12) enables the internal moment within the closed span to be related to the second 
derivative of the z-direction deformation (w) of the web as, 


<i 2 w(x) 

dx 2 


M(x) 


F T 

cs cs 


m d 



( 12 ) 


where E cs is the effective plane strain Young’s modulus of the closed span in the x-direction, and I cs is the 
area moment of inertia of the closed span face normal to the x-direction. Note that the plane strain 
Young’s modulus differs from the standard Young’s modulus by a factor of (1 - v xy v yx ) For a 
composite laminate, this value can easily be determined from classical lamination theory. Integrating 
equation (12) and imposing the zero rotation boundary condition at point D gives, 


<7w(x) 

dx 


F T 

cs cs 


M d x + P 



(13) 


which allows the rotation at point B (measured as positive in the clockwise direction) to be written as, 

dw(x = ) _ Q Mp hi | (14 

dx E cs I cs 6 E cs I cs 

The axial (x-direction) deformation at point B is due to the axial force in member BD, and can be written 
as, 


u B 


D . 


F A 

cs cs 


T _ C X~ F T 

Lcs E A L ° s 

CS I CS 


(15) 


where equation (2) has been employed. 

Next, equations for the rotation and x-direction displacement at point B will be developed by 
considering the deformation of members AC (crown) and AB (web). Equating these expressions with 
equations (14) and (15) will provide the two additional equations required to determine C x , M c , and M D . 

It is noted from figure 3 that, although the boundary condition at point C allows z-direction 
translation, no rotation is permitted at this point. Thus, the bending and rotation of member AC (crown) 
is due to the constant internal moment, M c , while the axial deformation is due to the axial force, C x (see 
fig. 7). Therefore, the rotation, z-direction displacement, and axial displacement at point A can be written 
as (ref. 12), 


q _ M c L cr 
A f T 

cr cr 


(16) 
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(17) 


w A = 


MAr 

IF T 

J ^ ± ^cr ± cr 


+ w c 


u A 


C x 

Ecr^-cr 


(18) 


where E cr is the effective plane strain Young’s modulus of the crown in the x-direction, and I cr is the area 
moment of inertia of the crown face normal to the x-direction, and w c is the z-direction displacement at 
point C. 

The web-crown combination (member BAC) is dismembered as shown in figure 8. A new, primed 
coordinate system is introduced for the web (member AB), with the x'-direction along the web (figs. 8 
and 9). The section of member AB, shown in figure 10, is used to determine the internal moment, M(x'), 
by applying equilibrium, 


= -> M c -M{x')+C x sin(e)x' = 0 (19) 

Solving for the internal moment, 


M(x') = M c + C x sin (0) x' (20) 

Beam theory (ref. 10) enables the internal moment within the web to be related to the second derivative of 
the z'-direction deformation ( w ') of the web as, 


z 



Figure 7. — Free body diagram of a section of 
member AC (crown). 



Figure 8. — Free body diagram of dismembered member BAC (crown and web). 
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X 


Q 


B 


B 
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Figure 9. — Rotated free body diagram of member AB 
(web). 

z' 


C x cos(0) 


M(x') 


C x sin(0) 

M C 


0 


a y c x cos(e) 

C x sin(0) 

Figure 1 0. — Free body diagram of a section of member AB. 


<i 2 w'(V) M(V) 

dx I w 


-i-[M c + C x sin(e)x'] 


( 21 ) 


where E w is the effective plane strain Young’s modulus of the web in the x'-direction, and I w is the area 
moment of inertia of the web face normal to the x'-direction. Integrating this expression provides the first 
derivative of w', which is equivalent to the rotation, 


dw’(x f ) 
dx ' 


F T 


M c x ' + ~C X sin(0)^ 


.'2 


+ C\ 


(22) 


where C\ is a constant of integration, which is determined by imposing continuity of rotations at point A. 
Noting that a positive internal moment in the web will result in a positive (clockwise) rotation, 


®A = 


dw f (x f = 0) 
dx' 


= Q 


(23) 


Equating equations (23) and (16), 


Q = 


M c^cr 

F T 

^ cr 1 cr 


Integrating equation (22) while employing equation (24) gives, 


w '( x ')=-¥£- x ' 2 + 


C x sin (e) x '2> + McLer x , + C2 


(24) 


(25) 


where C 2 is a constant of integration, which is found in terms of the web z' -direction deformation at point 

A, W A , 
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w'(x' = 0) = w' A =C 2 

Similarly, the axial deformation of the web is written as, 


(26) 


4*0 = %^ — x' + u' A 


E W A W 


(27) 


where is the web x'-direction deformation at point A. The rotation equations relating the web primed 
coordinate system deformation and the unprimed coordinate system deformations are given by, 


u — —u cos(0) + wsin(0) 
w = wsin(o)+ wcos(0) 


Therefore, from equations (25) to (28), we have 


w 


'h)= 


?F T 


— —x' 2 + McL cr + u ^ CO s(0) 


6E W I W 


F T 

^ cr 1 cr 


i(x) = CxC0S ( 9 ) x ' - Ua cos(0) + wj sin(0) 


Inverting equation (28) yields, 


F A 


u = -«'cos(0) + w'sin(0) 
w = «'sin(0)+ w'cos(0) 


(28) 


(29) 

(30) 


(31) 


which enables the x-direction deformation of the web at point B to be written as, 

ub = —u'b cos(0) + w'b sin(0) = -u '(. x ' = L w )cos(0) + w'(x' = L w )sin(0) 

= - — C ° S ^ L w +ua cos 2 (0) — wa sin(0)cos(0)+ (32) 

E W A W 'hEyyl yy 

+ c, sin 2 (9) ^ + mcLqa +UA Sin 2 (e)+WA cos( e) sin (e) 

6E W I W E cr I cr 

Recognizing that sin 2 (0) + cos 2 (0) = 1 and substituting for u A in equation (32) using equation (18) gives, 

(33) 


C x cos 2 (0) r , Mg sin(Q) r? , C x sin 2 (0) r , , M c 4 7 .sm(Q) r C x L cr 
B E A w 2 E I w 6E I w El w E A 

J -'XV JZL 'W j-^ rr ± rr ±^ rr ^L rr 


Equating equation (33) with the expression for u B determined from member BD, equation (15), 
C x cos 2 (0) r i M c sin(0) r2 , C x sin 2 (0) , , M c L cr sin(0) T C x L cr _C X -F T 


L L + -J^^±1± L 2 
FA 2 F J 


6 E W I w 


F T 

^ cr 1 cr 


LT 


F A 

^ cr ^ cr 


F T 

J ^cs ± cs 


(34) 


Noting that the internal moment within the web will cause a positive (clockwise) rotation at point B, 
equations (22) and (24) can be used to write the web rotation at point B as, 
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( 35 ) 


dw (x E w ^j M (jL w C x sin(0)A w M cL ( 

= Vb = ~ T + — r— + 


E W I W ^E W I W E cr I cr 


Equating equation (35) with the expression for 0 5 determined from member BD, equation (14), gives, 

McK , C * sm (9)4 | M c L tr M d L cs | PL] :s ( 

Tj 1 j 9/7 T F T 77 7 f\F T 

I - J w 1 w 1 w ^cr 1 cr rj cs 1 cs ® I - J cs 1 cs 

Equations (4), (36), and (34) now form three algebraic equations for the three unknown quantities, C x , 
M c , and M D _ Reiterating these equations, we have: 

Equilibrium 


M D+ M c + hC x + v[^-L cs \ + P\^L cs -^-\-M = 0 

Z y Z o 


Rotation Continuity at Point B 


McK , C x s m (0)Ll | M c L cr _ M D L cs | P I? cs 


F T 

w w 


2E W I w 


FT FT f>F T 

^cr^cr J ^cs ± cs KJJ ^cs ± cs 


x-Displacement Continuity at Point B 


Qco 7(e) l + Mcsinie ) L 2 1 C x sin 2 (e) ^ 3 | M c L cr sin (0) ^ C x Z cr V ~ ^ r 




2E W I w 


6E w J w 


F T 

cr cr 


w F A F A cs 

cr cr cs cs 


To solve these three equations, we rewrite equation (39) as, 


CJ- 


Z w cos 2 (0) I? w sin 2 (0) L c 


E w AL W 


6E W I w E cr A cr E cs A cs 


^sin(0) L w L cr sin(0) _ -FL CS 


2E W I w 


F T 

cr cr 


F A 

cs cs 


and identify the terms, 


K i=- 


L w cos 2 (9) | A sin 2 (9) 


FA fiF T FA FA 

UjC/ w 1 w ^ cr ^ cr ^ cs < cs 


L w sin(0) L w L cr sin(0) 

A 7 = 1 

IF T FT 

w w cr ± cr 


F A 

^ cs ^ cs 


Thus, from equation (40), we have, 
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(44) 


M c = 


A _r A 

k 2 x k 2 


Similarly, equation (38) can be written as, 


M r 


F J 

w w 


- + - 


F T 

^ cr cr 


K , 


+ c v 


A sin (0) 

2 E W I W 


= M 


D 


K a 


4. 


F T 

. ^cs ± cs 


+ - 


PLi 


^cs I CS 


Ks 


(45) 


with the following terms identified, 


a=- Zw 


A T 

w w 


- + - 


A T 

^ cr ± cr 


(46) 


ifa = 


^w sin ( e ) 

2 A 7 w 


£<=- Lcs 


F T 

^ cs ± cs 


(47) 

(48) 


P A 

*,=- cs 


67- tv f Ci 


(49) 


Thus, from equation (45) combined with equation (44) we have, 


A)=C, 


k 4 k x k 3 


k 5 k 2 k 5 


+ - 


*1^3 B 2 


k 2 k 5 k 5 


(50) 


Substituting equations (50) and (44) into the equilibrium equation (37) gives 

f c V 6 


C T 


k 4 k x k 3 


k 5 k 2 k 5 


+ - 


b,k 3 


. A + A _ c v Fl+hc +v 


k 2 k 5 k 5 k 2 


Ko 


A 

V 2 


+p 


s_ _s^ 

J c *~ A 


-M = 0 


(51) 


which can be solved for C x in terms of known quantities, 


M + P 


C v 


Vs 

v 8 2 c 


+ nf cs -7 
2/ 


A 

K, 


—2- + 1 
As . 


+ - 


A 


if. 


AA_A + /> 


if 5 if 2 if 5 


if. 


(52) 


The remaining unknowns, M c and M D , can then be determined from equations (44) and (50). 

The forces and moments in each component are given in terms of the three unknowns as follows, 


F = F -C 

1 cs 1 


V cs =Px 


(53) 

(54) 
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Px 

m cs =m d +— 

(55) 

V 

II 

(56) 

O 

II 

(57) 

M cr - M c 

(58) 

F W = C X cos(e) 

(59) 

F w = C x sin(0) 

(60) 

,=M c +C x sin(e) x' 

(61) 


where F t denotes an axial force and V t denotes a shear force. It is noted that the crown has constant axial 
force and moment and zero shear force, the web has constant axial and shear force and a linear moment, 
and the closed span has constant axial force, linear shear force, and quadratic moment. 

It should be noted that the areas and moments of inertia of each component are linearly dependent on 
the out-of-plane dimension of the hat (i.e., the dimension in the y-direction in fig. 3). However, examining 
equations (37) to (39), it is clear that this dimension drops out the governing equations. Thus, it can be 
concluded that if the applied moment, M, applied shear force, V, and applied axial force, F, are given as 
force and moment resultants rather than a moment and forces, and the distributed force due to pressure is 
given as a true pressure rather than a force per unit length, the force and moment solutions in equations 
(53) to (61) will be resultants as well. 

B. Symmetry Boundary Condition between Corrugations 

Next a symmetry boundary condition is considered at the midpoint between the panel corrugations, 
point H in figure 11. The corresponding free body diagram is identical to that shown in figure 4. In this 
case, the normal force and moment at point H are unknown reactions, while the shear force at point H is a 
known applied load. The primary distinction from the derivation in the previous section is that the 
deformation of the open span and face-sheet flange combination must now be considered in order to 
determine the unknown reactions at point H. 

A free body diagram of a section of member HGB is shown in figure 12, where j takes on the value of 
os or com depending on whether the section cut is taken within the open span or flange/facesheet 
combination. A new x coordinate direction has been introduced. Applying equilibrium, 

£f z =0 -> Vj = V - Px (62) 

?2 

Ya M Pt-a M j + P — -M + VjX = 0 (63) 


Substituting equation (62) into equation (63) gives, 
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Figure 1 1 . — Analysis geometry for a hat stiffened panel — symmetry at point H. 

z 
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s-tlTTTTTl^— ’ 

V P 

Figure 1 2. — Free body diagram of a section of member HGB. 


M,- =M + P Vx 

3 2 


(64) 


Once again utilizing beam theory (ref. 12), the second derivative of the z-direction deformation within the 
open span can be written as, 


d2w os (*) _ M os (*) 


1 


dx 


FT FT 

J ^os ± os J ^os ± os 


M + P Vx 

2 


(65) 


Integrating equation (65) and imposing the zero rotation boundary condition at point H ( x = 0 ) gives, 

dwps (*) 1 


dx 


F T 

^ os ± OS 


~3 ~2 

Mx + P- — V— 
6 2 


( 66 ) 


which allows the rotation at point G (measured as positive in the clock-wise direction) to be written as, 


dw os (x = L os ) 


ML„ 


dx 


- + - 


PLt 


VLt 


FT 6F T IF T 

^ os 1 os os 1 os os 1 os 


(67) 


Similarly, applying beam theory (ref 12) to the flange/facesheet combination gives, 
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d ^ com (-*0 _ M com (-*0 


dx 


F T 

com com 


F T 

com com 


M + P Vx 

2 


( 68 ) 


Integrating equation (68), 


d^com (-*0 _ Mx + Px 


Vx z 


dx 


FT f>F T IF T 

com com ^ com com com ± com 


-+c 


(69) 


where C is a constant of integration which can be determined by applying equation (69) to determine the 
rotation at point G (measured as positive in the clock-wise direction), 


d w com Aw ) _ p _ ML os 

— — — t) r7 — - 


PLt 


VLi 


dx 


- + ^ ^ + C 

FT f>F T IF T 

^ com com com com com com 


(70) 


and equating equations (70) and (67), 


C = 


ML„ 


- + - 


PLt 


VLt 


ML„ 


PLt 


- + - 


VLi 


F T ftF T 9 F T FT f\F T IF T 

^ os 1 os UjC/ os 1 os Arj os 1 os ^ com 1 com UjC/ com 1 com com 1 com 


(71) 


The rotation at point B (measured as positive in the clock-wise direction) can then be written as, 
d^com (^ — Aw + L com ) 


dx 


— — 6 R — 


3 2 

M (Aw + L com ) P^L os + L com ) V(^L os + L com ) 


F T 

com ± com 


- + - 


f>F T 

^ com ± com 


9 F T 

com * com 


+ C (72) 


Substituting equation (71) into equation (72), 


Q B =M 


FT FT 

os os com * com 


K. 


+ - 


( Aw + ^com ) Aw A 


T 3 

os 


F T 

com com 


F T 

^ os OS 


+ - 


(A os Aom ) Aw ^ VL ( 


F T 

com com 


F T 

^ os ± OS 


with the following terms identified, 


(73) 


* 6 =“ 


L n 


FT FT 

os ± os ^ com ± com 


(74) 


B 3 =- 


(L +L \ -I? F 

V os ^com ) os N 


F T 

com com 


F T 

os os 


+ - 


V 


(A os "I" Aom ) A is VL ( 


F T 

com com 


F T 

os os 


(75) 


Equation (73) can then be written as, 


e B = MK 6 +B 3 


(76) 
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while equation (14) can be written as, 


0 g -M d K 5 +b 2 


where equations (48) and (49) have been used. Equating equations (76) and (77) gives, 


M=M d 


K^ + B 2 
K* Kc 


h_ 

k 6 


(77) 


(78) 


Substituting equation (78) into equation (37) provides the replacement for the equilibrium equation in 
terms of the three unknowns, C x , M c , and M D , 


Equilibrium 


M 


D 


f v \ 

1_^5 

K 6 J 


+ M c + hC x + V 


— U 


+ P 


r s_ _ ^ 

J “"8, 


A + *L = 0 

K 6 k 6 


(79) 


Noting from figure 4 that the axial force in member HGB is constant and equal to F, the axial 
displacement at point G can be related to that at point B as, 


U G ~ U B ~ F 


F A 

V com com J 


(80) 


while the axial displacement at point H can be related to that at point G as, 


u h ~ u g~F 


L n 


V E os ^os J 


(81) 


Substituting the expression for u B obtained from the closed span, equation (15), into equation (80) and the 
result into equation (81) yields, 

u H = C x — ^ F-^ F—^om F _^ = 0 ( 82 ) 

FA FA FA FA 

^cs^^cs ^cs^^cs ^comcom ^os^^os 

where the zero axial displacement boundary condition at point H (see fig. 11) has been imposed. Solving 
equation (82) for F yields, 


F = C V 


jj cs 

F A 

^cs^cs 


cs | cum I us 

\ F A FA FA 

V cs^cs ^ com ^com ^os^os J 


K n 


(83) 


or 


f = c x k 7 


(84) 
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with 


Kh — 


F A 

cs cs 


L. 


L, 


_j cum |_ 

FA FA FA 

^cs^^cs ^ comcom ^os^^os 


Substituting equation (84) into equation (39) provides the replacement for equation (39), 


(85) 


x-Displacement Continuity at Point B 


C x c os 2 (e) z | M c sin (9) ^. 2 | C x sin 2 (e) z3 | M c L cr sin(e) z 


F A 


IF T 


6F T 


F T 

^ cr ± cr 


C I 
F A 

cr cr 


a 


F A 

^ cs ^ cs 


-L -F^l 

^ CS p A 

^ CS I cs 


( 86 ) 


The rotation continuity equation (38) remains valid. 

To solve equations (79), (86), and (38) for the three unknowns, C x , M c , and M D , equation (86) is 
rewritten as, 


C x K l+ M c K 2 = 


z£xKlLci 
F A 

cs cs 


(87) 


where equations (41) and (42) have been employed. Upon rearrangement, equation (87) becomes, 


C v 


K l + 


KjLcF 

F A 

^ cs i cs J 


+ M qK2 — 0 




K[ 


and replaces equation (44) when written in the form, 


m c = ^-cA 

k 2 K 2 


with, 


and, 


K\ = K, + KjLcs 
FA 

^ CS I CS 


B[ = 0 


Substituting equation (89) into equation (45), 

+ B[K^_B 2 
K 2 K 5 k 5 


m d =c x 


K 4 
K, 


K 2 K 5 


( 88 ) 


(89) 


(90) 


(91) 


(92) 


NASA/TM— 2008-2 1 543 8 


16 



Substituting equation (92) into equation (78), the unknown reaction moment at point H can be written as a 
function of C x , 


M = C\ 


K 4 K[K 2 
k 6 k 2 k 6 


+ 


B[K 3 

k 2 k 6 


B3_ 

k 6 


(93) 


Substituting equation (92), and (89) into the equilibrium equation (79), 


C v 


k 4 k[k 3 


k 5 k 2 k 5 


b[k 3 


+ - ! -^--— + — -C x — + hC x + V 

K 2 K 5 K 5 K 2 K 2 


+p 


— 4. 


-L cs 

s 2 ' 

-C r 

X 

K[K 3 " 

2 

V 

" 8 J 

X 

X 

*2*6 J 


b:k 3 b 3 n 

1 2 - H — = 0 


k 2 k 6 


K. 


(94) 


This equation replaces equation (51) and can be solved for C x . 


c x =- 


v 8 2 c 


+ V 


L 

^ CS 0 

2 j 


A 

X o 


*L + l3- 


V^5 


4 


+ - 


6j 


^2 g 3 
*5 




k 5 k 2 k 5 




x 


4, 


*2*6 


(95) 


Equation (95) replaces equation (52) for the present case where symmetry is specified at point H (see 
fig. 11) and the reaction normal force, F, and moment at point H, M, are unknown. Once C x is 
determined, F can be calculated from equation (84) and Mean be calculated from equation (93). 

The forces and moments in the open span and facesheet-flange combination are given by, 


77 _ 17 _ 77 

r com r os 1 


V -V 

r com r os 


V-Px 


(96) 

(97) 


M com M os M Vx + 


Px 


(98) 


It should be noted that equations (96) to (98) are valid for the applied force and moment boundary 
conditions described in section III.A. 


IV. Results and Discussion 

Five verification examples are presented below comparing the current solution within HyperSizer to 
NASTRAN finite element stress results. The first four cases consider hat stiffened and two sheet 
corrugated panels, with isotropic adherends, subjected to tensile and moment loads at point H. A final 
example then considers a composite hat stiffened panel subjected to internal pressure with the symmetry 
condition applied at point H. 

A. Example 1: Hat Panel Loaded with Tensile Force 

The first example problem considers a hat stiffened panel subjected to a tensile force resultant of 
1 lb/in transverse to the direction of the hat stiffeners. It is assumed that the hat stiffeners are bonded to 
the face sheet with an epoxy adhesive. The problem dimensions and isotropic material properties are 


NAS A/TM— 2008-2 1 543 8 


17 



summarized in figure 13. Figure 14 shows a plane strain NASTRAN finite element model of the problem, 
illustrating the deformed shape. The resulting shear stress, x vz , within the facesheet to flange joint region 
is shown in figure 15(a). Figure 15(b) compares the shear stress plotted along the middle of the adhesive 
(dashed line in fig. 15(a)) in the finite element solution with the HyperSizer Joints adhesive shear stress 
solution. As shown, the agreement is excellent. 



Unit weight (Ib/ft 2 ) 3.493 

Top face — Thickness (in.) 0.1 

Core web — Thickness (in.) 0.09 

Bottom face — Thickness (in.) 0 

Panel — Height (in.) 1.747 

Corrugation — Spacing (in.) 4 

Bottom crown — Width (in.) 0.8 

Core web — Angle (degrees) 82.24 


Top flange — Width, hat only (in.) 2.079 

Top clear span— Free Width (in.) 0.6969 

Top flange — Thickness (in.) 0.09 

Bottom crown — Thickness (in.) 

Facesheet and Hat: Aluminum 

E = 10 Msi 
G = 3.846 Msi 



To-scale hat stiffened panel 

Adhesive: Epoxy 

E = .445 Msi 
G = .165 Msi 
v = 0.35 

Thickness = 0.009 in. 


v = 0.3 


Figure 13. — Example 1 and 2 hat stiffened panel parameters. 



Figure 14. — NASTRAN finite element model deformed shape for example 1 hat stiffened panel. 


NAS A/TM— 2008-2 1 543 8 


18 






Figure 15. — Shear stress, x yz , (a) in NASTRAN finite element model (b) plotted along adhesive dashed line in 
(a)), comparing finite element model with HyperSizer Joints. 
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(a)l 



Figure 16. — Peel stress, a zz , (a) in NASTRAN finite element model (b) plotted along adhesive dashed 
line in (a)), comparing finite element model with HyperSizer Joints. 

Figure 16(a) shows the peel stress field in the finite element model solution, while figure 16(b) 
compares the peel stress along the adhesive (dashed line in fig. 16(a)) with the HyperSizer Joints solution. 
The agreement is excellent near the free edge of the hat flange, while in the web region, where the peel 
stress is lower, there is some discrepancy. This is expected as the HyperSizer Joints model considers only 
a bonded doubler (see fig. 1) and accounts for the presence of the stiffener web only as an effective 
boundary load. As the stiffener becomes thicker, and less shell-like, one would expect such effective 
boundary loads to be less representative of a continuum finite element solution. 

B. Example 2: Hat Panel Loaded with Applied Moment 

The second example problem considers the same hat stiffened panel, but now it is subjected to a 
moment resultant of-1 in.-lb/in. As before, the problem dimensions and material properties are 
summarized in figure 13. Figure 17 shows a plane strain NASTRAN finite element model of the problem, 
illustrating the deformed shape. The resulting shear stress, x yz , within the facesheet to flange joint region 
is shown in figure 18(a). Figure 18(b) compares the shear stress plotted along the middle of the adhesive 
in the finite element solution (as indicated by the dashed line) with the HyperSizer Joints adhesive shear 
stress solution. As shown, the agreement is excellent. 
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Figure 17. — NASTRAN finite element model deformed shape for example 2 hat stiffened panel. 



(a)l 



Figure 18. — Shear stress, x yz , (a) in NASTRAN finite element model (b) plotted along adhesive dashed line in (a)), 
comparing finite element model with HyperSizer Joints. 
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(b) - 250 J 
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Figure 19. — Peel stress, a zz , (a) in NASTRAN finite element model (b) plotted along adhesive dashed line in (a)), 
comparing finite element model with HyperSizer Joints. 

Figure 19(a) shows the peel stress field in the finite element model solution, while figure 19(b) 
compares the peel stress along the adhesive with the HyperSizer Joints solution. The agreement is 
excellent near the free edge of the hat flange, while, once again, in the web region, where the peel stress is 
lower, there is agreement is not as good. 

C. Example 3: Two Sheet Corrugated Panel Loaded with Tensile Force 

The third example problem considers a two sheet corrugated stiffened panel subjected to a force 
resultant of 1 lb/in. As shown in figure 20, a two sheet corrugated panel is similar to a hat stiffened panel, 
with the difference being the fact the two sheet has continuous flanges. The problem dimensions and 
material properties are summarized in figure 20. Figure 21 shows a plane strain NASTRAN finite element 
model of the problem, illustrating the deformed shape. The resulting shear stress, x yz , within the facesheet 
to flange joint region is shown in figure 22(a). Figure 22(b) compares the shear stress plotted along the 
middle of the adhesive in the finite element solution with the HyperSizer Joints adhesive shear stress 
solution. As shown, the agreement is excellent. 
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Figure 23(a) shows the peel stress field in the finite element model solution, while figure 23(b) 
compares the peel stress along the adhesive with the HyperSizer Joints solution. As shown, there is a 
significant discrepancy between the HyperSizer Joints solution and the finite element model in this case. 
While both solutions show a trough and a peak, the location of these features is not consistent. Again, the 
approximate way in which the presence of the web is accounted for in the HyperSizer Joints solution (as 
an effective boundary load) is the likely cause of this discrepancy. 



Unit weight (lb/ft 2 ) 3.749 

Top face — Thickness (in.) 0.1 

Core web — Thickness (in.) 0.09 

Bottom face — Thickness (in.) 0 

Panel — Height (in.) 1.747 

Corrugation — Spacing (in.) 4 

Bottom crown — Width (in.) 0.8 

Core web — Angle (degrees) 83.68 


Top flange — Width, hat only (in.) 2.855 

Top clear span — Free Width (in.) 0 

Top flange — Thickness (in.) 0-^9 

Bottom crown — Thickness (in.) 0.09 

Facesheet and hat: Aluminum 

E = 10 Msi 
G = 3.846 Msi 



Conceptual two-sheet stiffened geometry 


To-scale two-sheet stiffened panel 


Adhesive: Epoxy 

E = .445 Msi 
G = .165 Msi 


v = 0.3 v = 0.35 

Thickness = 0.009 in. 

Figure 20. — Example 3 and 4 two sheet corrugated stiffened panel parameters. 



Figure 21. — NASTRAN finite element model deformed shape for example 3 two sheet corrugated 
stiffened panel. 
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(b) -3 J 

y (in) 


Figure 22. — Shear stress, x yz , (a) in NASTRAN finite element model (b) plotted along adhesive (dashed line in (a)), 
comparing finite element model with HyperSizer Joints. 
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(a)l 



Figure 23. — Peel stress, a zz , (a) in NASTRAN finite element model (b) plotted along adhesive (dashed line in 
(a)), comparing finite element model with HyperSizer Joints. 


D. Example 4: Two Sheet Corrugated Panel Loaded with Applied Moment 

The fourth example problem considers the same two sheet corrugated stiffened panel, but now it is 
subjected to a moment resultant of-1 in.-lb/in. As before, the problem dimensions and material properties 
are summarized in figure 20. Figure 24 shows a plane strain NASTRAN finite element model of the 
problem, illustrating the deformed shape. The resulting shear stress, x yz , within the facesheet to flange 
joint region is shown in figure 25(a). Figure 25(b) compares the shear stress plotted along the middle of 
the adhesive (dashed line in fig. 25(a)) in the finite element solution with the HyperSizer Joints adhesive 
shear stress solution. As shown, the agreement is excellent. 

Figure 26(a) shows the peel stress field in the finite element model solution, while figure 26(b) 
compares the peel stress along the adhesive with the HyperSizer Joints solution. The agreement, while 
better than that exhibited by the two sheet panel subjected to the applied force resultant, still shows some 
discrepancy with the finite element results. 
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Figure 24. — NASTRAN finite element model deformed shape for example 3 two sheet 
corrugated stiffened panel. 




Figure 25. — Shear stress, x yz , (a) in NASTRAN finite element model (b) plotted along adhesive (dashed line in 
(a)), comparing finite element model with HyperSizer Joints. 
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(b) 

Figure 26. — Peel stress, o zz , (a) in NASTRAN finite element model (b) plotted along adhesive (dashed line in 
(a)), comparing finite element model with HyperSizer Joints 


E. Example 5: Composite Hat Stiffened Panel Loaded with Pressure 

The final example problem considers a composite hat stiffened panel subjected to a uniform pressure 
load of 10 psi. The problem dimensions and material properties are summarized in figure 27. The chosen 
composite is a typical graphite/epoxy fabric with all zero plies. This configuration was chosen for the 
example problem so that the same problem could be solved in a 2-D plane strain FEA. A composite with 
off-axis plies would require a 3-D finite element model for comparison, which was not constructed for 
this study. Figure 28 shows the plane strain NASTRAN finite element model for the problem. As shown 
in figure 29, a NASTRAN plate finite element model was also constructed for this problem to provide 
further verification. 

Figure 30 shows the plane strain finite element model solution to this problem, illustrating the 
deformed shape. For comparison, figure 31 shows the deformed shape resulting from the HyperSizer 
closed form method. The HyperSizer results for min and max deflection at the center of the open span 
(where the deflections are greatest) are within 2.5 percent of both FEA results. Figures 32 and 33 compare 
the moment resultant and shear force resultant solutions of the closed form HyperSizer method and the 
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plate finite element model. Clearly, the agreement is excellent. A comparison between the plane strain 
FEA and the current method for the adhesive shear stress, x yz , is shown in figure 34. Figure 35 shows a 
similar comparison of the peel stress in the adhesive. In both cases, the current method shows excellent 
agreement with the plane strain finite element model, even in the vicinity of the web. This better 
correlation of the peel stress in the web region (compared to the previous example results) may be due to 
the thinner stiffener used in this example. This thinner geometry tends to make the hat stiffener act more 
like a shell, which in turn makes the effective load boundary conditions used to account for the web more 
representative. 



Unit weight (lb/ft 2 ) 1.174 

Top face — thickness (in.) 0.098 

Core web — thickness (in.) 0.028 

Bottom face — thickness (in.) 0 

Panel — height (in.) 1.747 

Corrugation — spacing (in.) 4 

Bottom crown — width (in.) 0.8 

Core web — angle (degrees) 82.24 


Top flange — width, hat only (in.) 2.079 

Top clear span — free width (in.) 0 

Top flange — thickness (in.) 0.028 

Bottom crown — thickness (in.) 0.028 



Conceptual hat stiffened panel geometry 


To-scale composite hat stiffened panel 


Facesheet and hat: Graphite epoxy fabric 

Ei = 11 Msi; E 2 = 10.1 Ms 
G 12 = G 13 = G 23 = 0.76 Msi 
V 12 = 0.0363 

Ply thickness = 0.014 in. 

Facesheet layup 

7 Plies; [0 7 ] - total thickness = 0.098 in. 

Loading 

Pressure = 10 psi 


Adhesive: Epoxy 

E = .445 Msi 
G = .165 Msi 
v = 0.35 

Thickness = 0.009 in. 

Hat flange, web, crown layup 

2 Plies; [0 2 ] - Total Thickness = 0.028 in. 


Figure 27. — Example 5 composite fabric corrugated stiffened panel parameters. 
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Figure 28. — Plane strain FEA mesh for hat panel subjected to pressure. 



Figure 29. — Plate FEA mesh for hat panel subjected to pressure. 
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Figure 30. — Plane strain FEA model deformed shape for example 5 composite hat stiffened panel with applied 
pressure. 



Figure 31. — Deformed shape from the current closed form formulation (scale factor different than in fig. 28). 
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Moment resultant 
(in.-lb/in.) 




Figure 32. — Moment resultant (a) in NASTRAN plate finite element model (b) plotted along facesheet comparing the 
plate FEA with HyperSizer. 
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HyperSizer 


■ FEA Plate Model 



Distance along facesheet, in. 


Figure 33. — Shear force resultant (a) in NASTRAN plate finite element model (b) plotted along facesheet comparing 
the plate FEA with HyperSizer. 
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Figure 34. — Shear stress, x yz , plotted along adhesive, comparing finite element model with HyperSizer Joints. 




Figure 35. — Peel stress, o zz , plotted along adhesive, comparing finite element model with HyperSizer Joints. 
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V. Conclusion 


A method for the rapid analysis of facesheet to flange bonded joints in composite corrugated panels 
has been developed. The bonded joint analysis is based on the HyperSizer Joints methodology (refs. 4 to 
9), which is an analytical approach for determining the stress fields arising in the joint adherents and 
adhesive. In order to obtain the boundary loads specific to a given corrugated panel needed for application 
within the HyperSizer Joints analysis, a beam analogy model has been developed. This beam analogy 
model allows for a moment resultant and normal and shear force resultants, as well as pressure, to be 
applied, the magnitudes of which (for instance) can be obtained from a global finite element loads model. 
Thus, starting on this global scale, the loads can be localized and applied to the corrugated panel flange to 
facesheet bonded joint, and the local ply-level in-plane and intralaminar stress fields can be calculated. 
Ply-level margins of safety can then be determined based on ply-level failure criteria, and the viability of 
a given corrugated panel design can be evaluated. This new capability has been incorporated within the 
HyperSizer Structural Sizing Software (ref. 1), enabling a new level of fidelity in the software’s 
corrugated panel design capability. 

Results comparing the method to finite element analysis for monolithic and composite hat stiffened 
and two sheet corrugated panels were presented. It was shown that the method predicts joint adhesive peel 
stresses that match well with the finite element result in most cases. The predicted joint adhesive shear 
stresses match very well with the finite element results in all cases examined. The main discrepancy 
between the method and the finite element peel stress results occurred at the point where the web 
intersects the flange at an angle (see fig. 3). This angle intersection geometry is accurately captured in the 
finite element model, while in the joint analysis, an effective boundary load is applied at this point and the 
detailed local stress field induced by the web angle within the continuum finite element analysis is not 
captured. It is thus expected that the predicted stress fields would be more approximate in this region. 

This discrepancy was not observed in the final example presented, involving pressure applied to a hat 
stiffened panel with a thinner stiffener flange. In sum, the method has been shown to give quite accurate 
results in an efficient closed form suitable for structural optimization problems in which thousands of 
potential corrugated panel configurations must be considered rapidly. 
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